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In models of Pt 111 and Pt 100 surfaces in water, motions of molecules in the first hydration layer are spatially 
and temporally correlated. To interpret these collective motions, we apply quantitative measures of dynamic 
heterogeneity that are standard tools for considering glassy systems. Specifically, we carry out an analysis in 
terms of mobility fields and distributions of persistence times and exchange times. In so doing, we show that 
dynamics in these systems is facilitated by transient disorder in frustrated two-dimensional hydrogen bonding 
networks. The frustration is the result of unfavorable geometry imposed by strong metal-water bonding. The 
geometry depends upon the structure of the underlying metal surface. Dynamic heterogeneity of water on 
the Pt 111 surface is therefore qualitatively different than that for water on the Pt 100 surface. In both cases, 
statistics of this adlayer dynamic heterogeneity responds asymmetrically to applied voltage. 



The behavior of water adjacent to hydrated metal sur- 
faces is central to the physics underlying electrochem- 
istry, and as such there is much theoretical work devoted 
to this topic. References [THT41 are examples. Our re- 
cent work in this area^ highlighted the role of interfa- 
cial properties occurring over nanosecond timescales. In 
this paper, we examine this dynamics in further detail, 
considering models of Pt-water interfaces and showing 
structural relaxations in the water ad-layers proceed ex- 
ceedingly slowly via mechanisms that are spatially het- 
erogeneous. 

Strong metal- water bonding forces ad-layer waters into 
structures that are antithetical to favorable hydrogen 
bonding. Accordingly, the ad-layers contain defects. Mo- 
tions are one- to two-orders of magnitude faster in the 
proximity of these defects than elsewhere on the surface. 
Further, rearrangements of these defects require coor- 
dinated motions of several water molecules. This type 
of heterogeneous dynamics is common in glass-forming 
materials. Some of the quantitative methods used to 
successfully interpret suc h beha vior in molecular dynam- 
ics models of glass former J^^EIl used here to elucidate 
the nature of water dynamics at metal surfaces. 

Model 

As in our earlier work,^^^^ we carry out molecular dy- 
namics simulations with a surface model of Siepman and 
Sprik.^^ The model provides an empirical water- metal 
chemisorption potential. It includes effects of the metal's 
electronic polarizability and its response to applied volt- 
age across an electrochemical cell. For the model param- 
eters we consider, the electrodes most resemble platinum 
surfaceJ P^^ ' with water-metal binding energies around 
0.4eV and a metal lattice constant of 4. OA. Consequently 
when the electrode is exposed to a bulk liquid reservoir, 
water molecules adsorb to the metal atoms and adopt 



orientations that allow for hydrogen bond formation with 
molecules adsorbed to adjacent lattice sites. Snapshots 
of the ad-layers formed spontaneously from a slab of liq- 
uid water are shown in Panels (a) of Figs. [5] and [l] which 
refer to Pt 111 and Pt 100 surfaces, respectively. The 
snap shots show configurations of the ad-layer waters 
only. The liquid in contact with the ad-layer is not shown 
for purposes of clarity. Other pictures of this system can 
be found in Refs. 13 and 14, 

Depending on the surface geometry, adsorbed water 
molecules form qualitatively different hydrogen bonding 
patterns. For the 100 surface, the four-fold coordination 
and lattice spacing are commensurate with a variety of 
two-dimensional hydrogen bonding patterns. Surface wa- 
ter molecules adsorb to nearly every available surface site 
and vacant surface atoms are rare. The surface patterns 
that emerge are domains distinguished by relative orien- 
tations of the adsorbed water molecules. Different colors 
in Fig. [ija) help identify different domains. Transient 
disorder in the form of line defects between domains of 
dipole aligned molecules relax slowly, as we shall see. 

For the 111 surface, the six-fold lattice coordination is 
incommensurate with planar arrangements that would 
allow all molecules to form four hydrogen bonds. In 
this case, the particularly stable configurations of wa- 
ter molecules involve the formation of three hydrogen 
bonds where a given water molecule donates two and ac- 
cepts a single hydrogen bond in a triangular manner with 
molecules residing on alternating coordination sites. Blue 
coloring in Fig. [2|a) helps identify these patterns. These 
particular arrangements are facilitated by the presence 
of nearby surface vacancies that help eliminate hydrogen 
bond frustration and account for 15% of all surface sites 
on average. Since vacancy diffusion is slow, the spatial 
distribution of vacancies is heterogeneous over relatively 
long times, as shown below. 



Orientational mobility and fields 
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The chemisorption energy (i.e., the metal- water bond- 
ing energy) is sufficiently strong that lateral diffusion of 
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FIG. 1. An instantaneous configuration and dynamic heterogeneity of tfie water ad- layer on tfie Pt 100 surface. Tfie ad-layer 
is in equilibrium with adjacent bulk water (not shown), (a) Hydrogen bonding patterns showing heterogeneous distribution of 
line defects. Colors of underlying electrode atoms highlight distinct domains of ad-layer waters with specific hydrogen bond 
arrangements, (b) and (c) Instantaneous and time averaged mobility fields, g(a; t) and ^(a;tobs), respectively, with tobs = '?"p/3. 
The snap shots in (a) and (b) are taken at the midpoint of the trajectory that is averaged to produce (c). Color code for the 
mobility field is given by ^(a, t)Ax^ where Ax = O.lA. 



water is rare on picosecond timescales. As such, the dy- 
namics of an adsorbed water molecule are dominated by 
rotations. To characterize changes in molecular orienta- 
tions, we focus on the variables Ui(t), which denote the 
unit vector parallel to the dipole of the ith water molecule 
at time t. 

There are different contributions to the time depen- 
dence of Ui{t) that act on different time scales. In typi- 
cal trajectories, Ui{t) oscillates rapidly about a reference 
direction before jumping to new reference direction. The 
oscillations reflect small-amplitude vibrations and libra- 
tions of the water molecule, and the jumps (or inst anionic 
events) reflect changes in long lived arrangements of sev- 
eral molecules. The former generally have periods of 1 
ps or less, while distinct molecular arrangements gener- 
ally persist for more than 5 ps, and usually much longer. 
Indeed, as we detail later, the mean persistence time is 
20 ns and 1 ns for the water ad-layers on the Pt 100 and 
Pt 111 surfaces, respectively. 

As a result of this separation of time scales, it is useful 
to consider 
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The value of the coarse-graining time, Jt, should be large 
enough to remove most vibrational and librational con- 
tributions to Ui(t), and thereby highlight contributions 
from structural reorganization. For this purpose, we use 
St = 2 ps. Distribution functions graphed later in this 
paper show that this value is indeed suitable. 

In terms of these coarse-grained coordinates, an instan- 



taneous measure of mobility is given by 



(2) 



To the extent that it is non-zero, the configuration at 
time t-\- At differs from that at time t. To use this quan- 
tity to count relevant reorganization events. At must be 
both greater than the coarse-graining time, St^ and much 
less than the time to de-correlate enduring orientations. 
For this purpose, we use At — 10 ps. This choice is 
suitable for describing ad-layer mobility because 10 ps 
is one- to two-orders of magnitude shorter than the rele- 
vant mean persistence times. Spatial resolution of mobil- 
ity can be resolved with reference to a mobility field. Our 
specific choice mobility field focuses on ad-layer molecules 
and is coarse grained over a length ^. In particular, we 
define 

g(a, t) = ^g.(t)(/)(a - a,(t); 00(2^ - \z,(t) - |). (3) 

Here, a^(t) is the two-dimensional projection of the zth 
water-oxygen position onto the plane of the electrode, 
and Zi(t) is the projection perpendicular to the plane; 
is the mean of Zi(t) for water molecules in the ad-layer; 
Q{x) = 1 or for X ^ or X < 0, respectively; ^(a; ^) is a 
delta-like function broadened over length scale ^; we use 
a normalized Gaussian, 0(a; ^) ex exp(— a^/2^^)), and set 
i = 1.5A. 

Snapshots of ^(a, t) for the 111 and 100 surfaces are 
shown in Figs. [2|b) and[l](b), respectively. Both sur- 
faces display significant heterogeneous dynamics, mani- 
fested as regions of high mobility within a background 
of low mobility. For the 100 surface, the high mobil- 
ity regions have a definite directionality in that they are 
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FIG. 2. An instantaneous configuration and dynamic heterogeneity of tfie water ad- layer on tfie Pt 111 surface. The ad-layer 
is in equilibrium with adjacent bulk water (not shown), (a) Hydrogen bonding patterns: water molecules engaging in the 
preferred hydrogen bond pattern (see right-hand side of panel) have their underlying electrode atom colored blue, (b) and (c) 
Instantaneous and time averaged mobility fields, ^(a;t) and ^(a;tobs), respectively, with tobs = '?"p/3. The snap shots in (a) and 
(b) are taken at the midpoint of the trajectory that is averaged to produce (c). Color code for the mobility field is given by 
^(a,t)Ax^ where Ax = O.lA. 



aligned with the underlying lattice, and they are diffuse 
in that they encompass several neighboring molecules. 
Reorganization is therefore collective, occurring along ex- 
isting hydrogen bond chains. Further, by transiently dis- 
rupting the local hydrogen bond network, reorganization 
of one molecule facilitates reorganization of neighboring 
molecules. Such behavior is characteristic of glassy sys- 
tems in general.^ 

For the 111 surface, the high-mobility patches also have 
a directionality, though the underlying lattice dictates 
that motion does not propagate in a straight line but 
rather jigsaws back and forth along the lattice. Mobile 
regions are more localized on the 111 surface compared to 
the 100 surface, but motion is similarly facilitated as evi- 
denced by the connectivity of the mobile regions. Movies 
illustrating the time-evolution of dynamic heterogeneity 
for the water monolayers on the 100 and 111 surfaces 
can be viewed at http://youtu.be/cCYLPtck DWkl and 
http://youtu.be/dODDWuP-qTg, respectfully. 

The time scales over which immobile and mobile re- 
gions interconvert are the time scales of persistence. The 
mean persistence time, Tp, is the structural relaxation 
time.^ Over times smaller than Tp, patches of high mo- 
bility occur with higher frequency at particular regions of 
space. The persistent structural degrees of freedom con- 
tributing to these motions can be highlighted by time- 
averaging the spatial mobility field. 



q{a;tohs) = - — / dtq{a,t). 

^obs Jo 



(4) 



where tobs is an observation time. The time averaging 
causes mobile domains to become diffuse, as the higher 



frequency features disappear. Heterogeneity disappears 
completely as tohs/'^p becomes large. 

For the 100 surface, this time averaging elucidates the 
role of line defects on the surface as promoters of reorga- 
nization. Shown in Fig. [TJc) for tobs = '^p/S, integrated 
motion is primarily along lines that are fully connected 
and span the entire system. The linear structures run 
parallel to the surface close packing, and orthogonal to 
each other. The features can be readily identified with 
persisting line defects on the surface by analyzing the hy- 
drogen bonding patterns. Figure [TJa) shows the molec- 
ular configuration at the midpoint of the trajectory used 
to calculate Fig. TJc), where the underlying metal lattice 
has been colored for one of four hydrogen bonding pat- 
terns (also shown). The location of line defects formed at 
the boundaries between these different domains coincides 
with the regions high mobility. 

For the 111 surface, time averaging the mobility field 
yields similar results. As shown in Fig. [2|c) for tobs = 
Tp/S, integrated motion exhibits domains of high mobil- 
ity. These domains, like the ones on the 100 surface, 
become more connected although their structure is not 
as simple. The locations of the domains are not directly 
related to the location of surface vacancies but rather are 
spatially anti-correlated to the presence of stable hydro- 
gen bond configurations. Molecules with this preferred 
hydrogen bonding pattern typically reside within inactive 
patches. Figure [2|a) is taken from the midpoint of the 
trajectory used for Fig. |2|c) and designates these partic- 
ularly stable hydrogen bond configurations as blue lattice 
sites. 
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Correlation and distribution functions 



function, 



These heterogeneous, facihtated dynamics described 
above dictate that the timescales associated with relax- 
ing a given tagged molecule will be large relative to the 
bulk liquid, and strongly dependent on its environment. 
The average timescale governing orientational reorgani- 
zation can be extracted from the dipole autocorrelation 



CW = {u(0)-u(t))-Ku)|2, 



(5) 
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FIG. 3. (a) The dipole autocorrelation function for water 
molecules adsorbed to the 100 (dotted red line) and 111 (solid 
blue line) electrode. The corresponding quantity for molecules 
in the bulk liquid is plotted as a dashed black line, (b) The 
probability distributions for the persistence, tp, and exchange 
times, tx, for the 111 surface, shown on a log scale where 
P(log t) = tP{t). (c) The probability distributions for the 
persistence and exchange times for the 100 surface. 



where the angle brackets represent an average over all 
electrode-adsorbed water molecules. Note that while (u) 
vanishes in the bulk liquid, it does not necessarily at the 
surface. As shown in Fig. [sj^a), water on both surfaces 
undergo orientational relaxation over timescales that are 
of 2 to 20 ns. The more ordered 100 surface yields signif- 
icantly slower relaxation than for that of the locally over- 
coordinated 111 surface. These times are 3-4 orders of 
magnitude larger than the characteristic relaxation time 
for bulk water, Tbuik = 5 ps. (For the computation of 
C{t) for bulk water, the averaging implied by the angle 
brackets is carried out over all non- ad-layer waters.) 

The functional form of the decay of C{t) for bulk wa- 
ter is an expected exponential decay, indicative of an 
uncorrelated Poisson process characterized by a single 
timescale, Tbuik- The functional form for the decay of 
C{t) for the adsorbed water, however, is more compli- 
cated and resembles the stretched exponential decay as- 
sociated with correlation functions of glass-forming liquid 
systems. Such deviations from Poisson statistics high- 
light the important temporal correlations within the ori- 
entational dynamics. 

These temporal correlations are manifested explicitly 
in the distributions of persistence and exchange times 
Persistence times, tp, are the waiting times for reorien- 
tations, and exchange times, tx, are the times between 
such motions. In cases of uncorrelated dynamics, persis- 
tence and exchange times are drawn from the same dis- 
tribution, fully characterized by the decay time of C{t). 
However, if the dynamics are temporally correlated then 
the probability for undergoing motion is conditionally de- 
pendent on what has already occurred, and these distri- 
butions will be therefore different.'^ In glassy dynamics, 
where motion in a region of space facilitates further mo- 
tion in neighboring regions of space and time, the most 
probable is much smaller than the most probable tp. 

This decoupling of exchange and persistence is clearly 
evident in Figures [sjb) and (c), which depict distribu- 
tions of persistence and exchange times for both the 111 
(b) and 100 (c) surfaces. The times are calculated by 
recording the time for a dipole to change its orientation 
by 50^ either from an arbitrary time origin, for the per- 
sistence time, or given a reorientation just occurred, for 
the exchange time. This angular displacement is chosen 
to be much greater than the librational motion of wa- 
ter on the surface and commensurate with reorientations 
between stable hydrogen bonding patters on the surface, 
as dictated by the lattice geometry. In order to avoid 
counting reorganizations that only occur transiently, we 
also stipulate that the rotation persists for an additional 
At. (For reference to similar calculations of exchange and 
persistence distribution, but for a model of a structural 
glass- forming liquid, see Ref. 23 ) 

For ad-layers to both metal surfaces, the mean ex- 
change time, Tx, is shorter by over an order of magnitude 



5 



from the mean persistence time. The mean persistence 
time, Tp, agrees with the 1/e time from the decay of C{t). 
This separation of timescales, reflecting a dynamic facil- 
itation mechanism, occurs in here because transient dis- 
order within the adsorbed monolayer gives rise to spatial 
variations in the relative hydrogen bond stability, which 
in turn facilitates heterogeneous dynamics. 

The total orientational mobility within the adsorbed 
monolayer is given by 

1 ^ 

where the sum is carried out over all N electrode- 
adsorbed molecules. {N can fluctuate due to exchange 
between the bulk and ad-layer, but the resulting time 
dependence of N has a negligible effect on Q{t) for the 
conditions we have studied.) This total orientational mo- 
bility variable serves as an order parameter distinguishing 
ad-layers of different mobility. Figure [4] shows the dis- 
tribution functions of of this order parameter, p{Q) for 
the 100 and 111 electrode surfaces. Both distributions 
are Gaussian around the mean, but both also possess fat 
tails at larger values of Q. The non-Gaussian tails in- 
dicate that highly active configurations are much more 
probable than one would expect given Gaussian statis- 
tics. The tails arise through the correlated nature of the 
surface dynamics. Their presence indicates, in principal, 
that it is possible to drive the ad-layers out of equilib- 
rium into distinct phases, one of high mobility and one 
of low mobility. 

This behavior does not change qualitativey under ap- 
plied voltage. Interestingly, the distributions of total mo- 
bility respond asymmetrically at the positive and nega- 
tive electrode. Figure |4] shows p{Q) for both electrode 
geometries and for the positive (Vq = 1.36V) and neg- 
ative {Vo = —1.36V) constant potential electrode. For 
both surfaces, the positive values of the applied poten- 
tial have the effect of shifting the distribution of mobility 
to lower values of Q. This shift to lower Q manifests as 
a slowing in the timescales for relaxation dynamics. At 
the negative potential electrode we observe the opposite 
effect, a shift in p{Q) to larger values of Q, indicating an 
increase in the surface relaxation times. 

The asymmetric behavior can be understood in light 
of our previously published work in which we demon- 
strated that for moderate values of an applied electrode 
potential the orientations of adsorbed water molecules re- 
spond asymmetrically with respect to applied electrode 
potential. The asymmetry arises because at the nega- 
tive electrode adsorbed water molecules can adopt an al- 
ternative orientation in which one oxygen- hydrogen bond 
is pointed straight towards the electrode thereby gaining 
a favorable image charge interaction between the partial 
positive charge on the hydrogen and the electrode. There 
is no corresponding orientation at the positive electrode 
and so at the negative electrode one relaxation pathway 
(rotation of an oxygen- hydrogen into the electrode) is 
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FIG. 4. (a) The probability distribution for the total ori- 
entational mobility, Q, plotted for the 100 electrode surface 
(dashed red line) and the 111 electrode surface (solid blue 
line), (b-c) The probability distribution p{Q) at different 
values of the applied electrode potential. Dotted black line 
correspond to the results at zero applied potential. Red and 
blue lines correspond to the results at the negative electrode 
(Vb = —1.36V) and positive electrode (Vb = 1.36V) respec- 
tively. 



enhanced while at the positive electrode the same path- 
way is inhibited. This asymmetry was found to be more 
prominent for 111 versus 100 due to the relative instabil- 
ity of the hydrogen bond network over the 111. 

Consistent with the asymmetric behavior predicted 
from our model, others have noted that the electronic 
structures of metal interfaces respond sensit ively to the 
electrostatic fields of the surrounding solvent .^^^ In par- 
ticular, electronic structure theor}^^^ indicates that dipole 
reorientations of adsorbed water molecules on a Ft 111 
surface can result in shifts in the potential of zero charge 
of up to 3 eV. The heterogeneous relaxation elucidated 
herein are therefore expected to modulate catalytic ac- 
tivities of metal surfaces. For a platinum-like surface, we 
expect modulating domains that extend over nanometers 
and that reorganize tens of nanoseconds. It remains to 
be seen whether asymmetric response can be tapped as 
a means to drive electrodes out of equilibrium in a fash- 
ion that will produce a non-equilibrium transition be- 
tween phases of high and low orientational mobility, and 
whether such a transition can be of practical use. 
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